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ABSTRACT 

We interpret the relation between the polarizing efficiency Pniax/E{B — V) and the wave¬ 
length of the maximum polarization Amax observed for 17 objects (including 243 stars) sepa¬ 
rated into two groups: “dark clouds” and “open clusters”. The objects are assigned to one of 
the groups according to the distribution of the parameter Amax- We use the model of homoge¬ 
neous silicate and carbonaceous spheroidal particles with the imperfect alignment and a time- 
evolving size distribution. The polarization is assumed to be mainly produced by large silicate 
particles with the sizes ry ^ rv,cut- The models with the initial size distribution reproducing 
the average curve of the interstellar extinction fail to explain the values of Amax ^ 0.65 /rm 
observed for several dark clouds. We assume that the grain size distribution is modified due 
to accretion and coagulation, according to the model of Hirashita & Voshchinnikov (2014). 
After including the evolutionary effects, Amax shifts to longer wavelengths on time-scales 
^ 20(nH/10^cm“^)“^ Myr where nn is the hydrogen density in molecular clouds where 
dust processing occurs. The ratio P„iax/E{B — V) goes down dramatically when the size 
of polarizing grains grows. The variations of the degree and direction of particle orientation 
influence this ratio only moderately. We have also found that the aspect ratio of prolate grains 
does not affect significantly the polarizing efficiency. For oblate particles, the shape effect is 
stronger but in most cases the polarization curves produced are too narrow in comparison with 
the observed ones. 

Key words: polarization — dust, extinction — galaxies; evolution — galaxies: ISM — ISM; 
clouds 


1 INTRODUCTION 

Non-spherical dust grains are responsible for the interstellar polar¬ 
ization phenomenon. The polarization becomes significant when 
interstellar grains are well aligned and have sizes comparable to 
the wavelength of the incident radiation. Alignment of dust parti¬ 
cles can arise due to their specific magnetic properties giving them 
the ability to efficiently interact with the interstellar magnetic fields. 

The interstellar linear polarization is characterized by the po¬ 
larization degree P and the position angle 9e or Oq measured either 
in the equatorial or galactic coordinate system, respectively. The 
polarization degree usually tends to a maximum in the visual and 
gradually decreases to the ultraviolet and infrared. The wavelength 
dependence of pola rization P(A) is d escribed by an empirical for¬ 
mula suggested in ISerkows^ ( Il973h and now called Serkowski 
curve 

P(A)/Pmax = exp[-Aln"(Amax/A)]. (1) 

* E-mail: n, voshchinnikov@spbu.ru 


This formula has three parameters: Pmax is the maximum degree 
of polarization, Amax the wavelength corresponding to it, and K 
the coefficient characterizing the width of the Serkowski curve. 
The values of Pmax in the diffuse interstellar medium usually 
do not exceed 10%, a nd the average value of Amax is 0.55 pm 
dSerkowski et aljri975h . The parameter K is related to the half¬ 
width of the normalized polarization curve 

W = Amax/A— — Amax/A+, (2) 

where A_,A+ are such that P(A+) = P(A_) = Pmax/2 and 
A_ < Amax < A+. The relation between W and K is as follows: 

W = exp[(ln2/A')^/^] - exp[-(ln2/A)^/^]. (3) 

Initial ly, the parame t er K w as chosen to be 1.15 dSerkowskil 1973lj . 
Later, IWhittet et alj (Il992h found a dependence of K on Amax in 
the Milky Way 

K = (1.66 ± 0.09)Amax -f (0.01 ± 0.05), (4) 

where Amax is in microns. 

The ratio of Pmax to the colour excess E{B — V) or visual 
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Table 1. Tai'get objects. 


N 

(1) 

Object 

(2) 

-A^stars 

(3) 

1 

(4) 

b 

(5) 

D, pc 

(6) 

(Aniax)i /J-m 

(7) 


do,a ± o-gg Q 

(9) 

References 

(10) 

\e(B-V)) 

(8) 






Dark clouds 




1 

Taunas: cloud 1 

31 

174 

-14 

147 

0.570 ± 0.049 

6.01 ± 2.09 

156.5 ± 13.1 

[1,2,3] 

2 

Chamaeleon I 

25 

297 

-15 

196 

0.594 ± 0.064 

9.46 ± 2.53 

101.3 ± 10.8 

[4,5,6] 

3 

p Oph: cloud 1 

10 

353 

+18 

133 

0.684 ±0.110 

3.85 ± 1.56 

31.4 ± 14.4 

[7,8,9,10] 

4 

cloud 2 

20 




0.659 ± 0.059 

6.19 ± 3.21 

98.3 ± 11.2 


5 

R CrA: cloud Al 

9 

359 

-18 

148 

0.758 ±0.104 

3.48 ± 1.62 

132.0 ±24.1 

[4,5,11] 

6 

cloud A2 

5 




0.734 ± 0.080 

2.21 ± 1.12 

52.1 ± 14.3 


7 

cloud B 

5 




0.690 ±0.071 

5.48 ± 2.57 

153.4 ± 11.0 


8 

CMa Rl: cloud 1 

5 

225 

-1 

690 

0.637 ± 0.030 

4.80 ± 2.06 

132.4 ± 13.0 

[12] 

9 

cloud 2 

6 




0.593 ±0.149 

2.31 ± 1.06 

75.7 ± 24.3 






Open clusters (and some dark clouds) 



10 

Taurus: cloud 2 

13 

174 

-14 

147 

0.556 ± 0.043 

5.76 ± 1.61 

26.1 ± 20.4 

[2,3] 

11 

a. Per 

36 

147 

-6 

213 

0.543 ±0.051 

8.10 ±2.63 

86.7 ± 13.8 

[13] 

12 

Musca 

16 

300 

-9 

171 

0.577 ± 0.017 

12.62 ± 2.98 

102.9 ± 3.9 

[4] 

13 

Southern Coalsack 

14 

301 

-1 

174 

0.570 ± 0.040 

6.31 ± 2.73 

75.3 ± 17.4 

[4] 

14 

NGC 654 

6 

129 

-0 

2410 

0.555 ± 0.019 

4.14 ±0.71 

87.3 ±4.5 

[14] 

15 

IC 1805 

16 

134 

+1 

2400 

0.540 ±0.031 

5.70 ± 1.14 

96.1 ±5.0 

[15] 

16 

NGC 6124 

11 

341 

+6 

563 

0.574 ± 0.039 

3.40 ± 0.84 

54.0 ± 4.7 

[16] 

17 

Cyg OB2 

15 

80 

+1 

1700 

0.473 ± 0.073 

2.07 ± 0.71 

146.6 ± 26.2 

[5,7] 


References: [1] -^firaov2009; [2] -^hittetetaL_200^; [3] -Ifst^222^; -^ndersson^Potter^OOV; [5] -^hittetet^^J992; 
[6] - McGregor et al. 1994: [7] - Martin et al. 1992: [8] - Wilkinn et al. 1982: [9] - Vrbaetal. 1993: [10] - Snow et al. 2008: 

[11] - Vrbaetal. 1981: [12] - Vrba et al. 1987: [13] - Coyne et al. 1979: [14] - Medhi et al. 2008: 

[15] -^edlue^af_2007; [16] -^rgneet^;_20W. 


extinction Av is called the polarizing efficiency. There exists an 
empirical limit to this ratio dSerkowski et al.ll975ll 


line of sight. As a basic model we use that of imperfectly aligned 
spheroidal particles in a regular magnetic field that was earlier ap- 


E{B - V) 


^ 9%/mag. 


(5) 

larization curves in a wide spectral range i 

Voshchinnikov & Dad 


I 2 OO 8 I: iDas et al.l 1201(1 ISiebenmoreen et al 

2014|). We consider 


Av RvE{B - V) 


^ 3%/mag. . 


( 6 ) 


where Rv is the total to selective extinction ratio. 

A qualitative explanation of the relation between the 
width of the interstellar polarization curve and the position of 
its maximum is connected to the fact that dust grains grow in 
the accretion and coagulation processes, which le ads to a nar¬ 
rowin g of the particle size distribution (see, e.g., IWhittet et al.l 
1 19921) . A c|U3.ntitsti vc int6rpr6t3.tion of the d epend ence K {Amax) 
was suggested by lAannestad & Gree nberd (1 19831) w ho co nsid¬ 
ered cylindrical particles and by Voshchinnikov et ^ ( 1201 3h and 
IVoshchinnikov & Hirashital ^201^ who used spheroidal grains. 

However, all three parameters of the Serkowski curve have 
been determined for a limited n umber of stars (less than 200; 
IVoshchinnikov & Hirashital l2014h . This is because one needs to 

S ierform observatio ns in more than four or even five bands 
Whittet et al.l 19921) to find the width of the polarization curve (pa¬ 
rameter K) with a good accuracy. In contrast, the parameters Pmax 
and Amax can be determined with a sufficient accuracy when stars 
are observed just in three or four bands. Hence, there exists a sig¬ 
nificant number of lines of sight for which both the polarizing ef¬ 
ficiency and the wavelength of maximum polarization have been 
estimated. 

The polarizing efficiency in a given direction depends on the 
size of polarizing and non-polarizing grains, the particle shape 
and the degree and direction of grain alignment. In general, the 
dust grain characteristics and the magnetic field can vary along the 


the time evolution of the grain size distribution due to accr etion 
and coagulation processes dHirashita & Voshchinnikqt ]|201'^ and 


involved new op ti cal constants of grain mate r ials l fonej 20121 ; 
Ijones et ^ 1201^ . IVoshchinnikov & Hirashital ( 12014 hereafter 
VH14) used such a model to investigate the relation between K 
and Amax that are the parameters of the Serkowski curve which 
characterize its width and maximum position, respectively. 

In this paper, we focus on interpretation of the observed de¬ 
pendencies of P-max/E{B — V) on Amax and show that the most 
important factor influencing these dependencies is the size of po¬ 
larizing grains. The paper is organized as follows: Sectionsl^and 
[^give a description of the observational data and the model used. 
Section |4] presents results of our modelling of the polarizing effi¬ 
ciency and their discussion, and Sectionj^contains the conclusions 
made. 


2 OBSERVATIONAL DATA 

We have collected available observational data on the polarizing ef¬ 
ficiency Pmax/E{B — V) and the maximum polarization position 
Amax from the literature. We started with a list of 160 stars with 
the known values of the parameter K considered in VH14. From 
the list, we selected stars closely located in the sky and with sim¬ 
ilar values of the polarization position angle. We extended these 
groups by including stars with the known values of Pmax, Amax, 
and E{B — V) and excluded the stars with significant wavelength 
rotation of the position angle and large errors in observational data. 
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Figure 1. Polarizing efficiency Pmax/E{B — V) versus the wavelength of maximum polarization Amax for dark cloud 1 in Taurus (a) and cluster a Per 
(c). The average values ^ e^b—^) '} plotted with standard deviations. Vertical and horizontal dashed lines correspond to the mean interstellar value of 
Amax = 0.55 fj,m and the upper limit of the polaiizing efficiency given by equation Panels (b) and (d) show the histograms of the distributions of Amax 
and E{B — V) for Taurus cloud 1 and a Per cluster, respectively. 


Morever, we chose only the probable members of the stellar clus¬ 
ters. 

Our final list contains 243 stars associated with 17 objects that 
are either a dark cloud or a stellar cluster. We separate the objects 
into two groups according to the patterns seen on their dependence 
of Pin 3 .x/E{B — V) on Amax (see Fig.[T}. All the objects are listed 
in Table [T] that includes the number of stars considered, the co¬ 
ordinates and distances to the objects, the average values of the 
wavelength of maximum polarization with the deviations. The last 
column gives the sources o f the observati onal data. We took the 
distances to the objects from iKnu d3 j2010ll and the papers cited in 


Table[T] Most of the objects are not very distant and either belong 
to the Gould Belt or are close to it. The exceptions are NGC 654 
and IC 1805 located in the direction of the anti-center of the Galaxy 
where the interstellar medium is transparent up to large distances 
jLallement et alj2014t) . 

For the clusters, polarization of their stars is mainly caused by 
dust in the diffuse medium in the line of sight. On the plots we see 
a cloud of points around some mean polarization efficiency, with 
Amax being usually in the interval 0.5 - 0.6 /rm (see Fig. ,d). The 
clusters produce what we call the open cluster pattern. 

For the dark clouds, when a part of the stars with measured 
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Figure 2. Number distribution of the position angles around the mean value 
(9 o,g fot ihe dark cloud 1 in Taurus and cluster a Per. Curves show the 
Gaussian fits of the observational data. 


polarization are seen through outer layers of the cloud, we meet 
another pattern called the dark cloud pattern. The stars observed 
mainly through diffuse medium, i.e. foreground, aside, distant stars, 
produce again nearly the open cluster pattern, but the stars seen 
through the dark cloud sometimes add a trend: the larger Amax, the 
smaller the polarization efficiency (see Fig. [Ili,b). 

Obviously, when the part of stars seen through the dark cloud 
is very small, we approach the open cluster pattern, even though the 
stars are at different distances, while for the clusters the distances 
are nearly the same. This pattern is seen in our data for Taurus cloud 
2, Musca cloud, and Southern Coalsack. 

The way we associate the stars with the objects is to some 
extent arbitrary and is based on the dat a published in the papers 
listed in Table 1 (e.g.. IVrba et alj[l9^ set off the cloud and 
cloud B in the dark cloud R CrA). In Taurus, two groups of stars 
with a relatively uniform distribution of the position angles (cloud 
1, (Iq = 145° - 175° and clou d 2, Oq = 2° - 40°) were de¬ 
scribed ^shchinnikovj2oI3) on the basis of the initial anal¬ 
ysis of [Messenger et alJ (Il997h . Stars in “clouds” are distributed 
around complexes of molecular clouds and partly are projected on 
them. In these cases, the interstellar extinction can differ signifi¬ 
cantly, in contrast to stars of “clusters” when it lies in a rather nar¬ 
row range (cf. Figs. [Ih and mi A negative correlation between 
Praa.^/E{B — V) and Amax for the dark cloud 1 in Taurus and a 
uniform distribution of data for a Per cluster are well seen in the 
figures. Note that neither such a correlation nor its absence is typ¬ 
ical of the clouds and clusters. Probably, the correlation is formed 
by local conditions and is partly a result of observational selection. 

Column (8) in Table [T] shows the average values of the 
polarizing efficiency derived from observations. These val¬ 
ues do not take into account the depolarization effect due to 
spatial variations of the magnetic field B. This effect can be 


^ Besides this, we distinguished clouds A1 and A2. 



Figure 3. Average polarizing efficiency {Pma-x./— V)) versus the 
average wavelength of maximum polarization (Amax) for dark clouds 
(squares) and open clusters (triangles) in the Milky Way. Vertical and hor¬ 
izontal dashed lines con'espond to the mean interstellar value of Amax = 
0.55 /rm and the upper limit of the polarizing efficiency given by equa¬ 
tion The number of stars, coordinates of objects, sources of data, etc. 
are presented in Table[T] 


found as the dispersion in directions of polarization angle. 
A similar approach has been applied to optical polarization 
jMvers & Goodmaj Il99ll : Ijones et ^ Il992h and submm polar¬ 
izatio n measur ed by Planck l Planck Collabo ration Int. XXXl j 
20161: _ I Planck Collaboration Int. l2016l: 


Planc k Collaboration Int. X Liyll2016l:lGuillell2015h . 

iMvers & Goodmanni99l[) studied the spatial patterns and dis¬ 
tributions of the direction of the interstellar polarization for several 
clouds, clusters and complexes of clouds. They found that the num¬ 
ber distributions often had a single local maximum with the disper¬ 
sion of 0.2 - 0.5 radians and the clouds had a larger dispersion 
than the clusters. We calculate the mean value of the position an¬ 
gle 6 ^ 0,0 and its dispersion agg for our objects (see column 9 in 
Table [Til. Figure [^presents the number distribution of the position 
angle for the dark cloud 1 in Taurus and cluster a Per together with 
the Gaussian fits. Note that the dispersions ^ are very close to 
the widths of Gaussian curves (13°1 vs 12°9 and 13°8 vs 13°6 for 
dark cloud 1 in Taurus and cluster a P er, respectively). Though w e 
used the data different from those of iMvers & GoodmanI (1 199 ll) . 
the dispersions found for a Per and Cha I do not differ more than 
by 20%. For the regions in Ophiuchus and Taurus the difference is 
about 50% mainly becaus e we considered separate clouds, while 
IMvers & GoodmanI ( Il99ll) large complexes. 

Figure |3 shows the average polarizing efficiency versus the 
average wavelength of maximum polarization for target objects. 
One clearly sees two patterns in distribution of the (Amax) values: 
for the open clusters and some clouds, the mean values of Amax 
are grouped around the average interstellar value equal to 0.55 ^m 
with rather small deviations, while for most dark clouds, the val¬ 
ues of (Amax) may be significantly larger than 0.55 pim and have a 
wider scattering. Note also that in two cases (Cha I and Musca) the 
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average polarizing efficiency exceeds the upper observational limit 
given by equation 0- 

We singled out the association Cyg OB2 which still remains 
one of the most actively s tudied objects in the Galaxy (see a his¬ 
tory of its investigations in ichentsov et alj2013ll . The polarimetric 
studies of stars in Cyg OB2 have demonstrated a low polarizing ef¬ 
ficiency and unusually small values of the maximum polarization 
wavelength. Therefore, Cyg OB2, marked by a cross in Fig. [3] oc¬ 
cupies an isolate position in the lower left corner. It is well known 
that in this direction the line of sight is aligned along the spiral 
arm, i.e. we should expect rather small polarization because of a 
moderate inclination of the magnetic field. From other side, the as¬ 
sociation Cyg OB2 is rather distant. So, depolarization of stellar 
radiation may occur due to quite different polarization in several 
clouds in the line of sight. Taking into account a contribution of 
foreground pol arization seems to lead to an increase of the polariz¬ 
ing efficiency jMcMillan & TapiJl977l : IWhittelll2015h . Evidently, 
polarization in Cyg OB2 direction needs a further analysis and dis¬ 
cussion. 

Undoubtedly, it would be interesting to compare the opti¬ 
cal polarization of our objects with the submm polarization mea¬ 
sured by Planck. Unfor tunately, the published results were com- 
puted at 1° resolution jPlanck Collaboration Int. XIXll2015h . So, 
the data can only suggest that this polarization is larger in the region 
Chamaeleon-Musca than that in Ophiuchus, Taurus and R CrA. 


3 MODEL 


Our model of interstellar dust grains is exactly the same as used by 
VH140We consider a mixture of silicate and carbonaceous homo¬ 
geneous spheroids with certain distributions over sizes and orienta¬ 
tions, but a fixed ratio a/b, where a and b are the major and minor 
spheroid semiaxes, respectively. Both the silicate and carbonaceous 
particles contribute to extinction, but the interstellar polarization is 
assumed to be produced mainly by silicate particles. This assump¬ 
tion is supported, in particular, by a correlation between the ob¬ 
served interstellar polarization degree and the abundance o f silico n 
in dust grains found in the work of IVoshchinnikov et al.l ilOl'j) . 
We also suggest that small silicate particles are randomly oriented 
when their sizes rv < rv.cut, where rv is the radius of a sphere 
whose volume is equal to that of the spheroid, rv = \fat? for 
prolate spheroids and rv = \faFb for oblate ones. 

Let us consider a dust cloud with the uniform magnetic field. 
The angle between the line of sight and the magnetic field is de¬ 
noted by fl (0° 5(12^: 90°). The extinction and linear polariza¬ 
tion of unpolarized stellar radiation produced by aligned rotating 
spheroidal particles are 


D '^V,max,j 


A{\) = 1.086 


E/ / 


ajlbj,rv, A, Q.) 


0 


xnj{rv) drv dl, (7) 


P(A) = 


D ’’V,max,Si 

C'pol,Si(^A,Si, aSi/^Si, , A, H) 

0 T'v,cut,Si 



xnsi{rv) dry dl x 100%, 


(8) 


^ In previous modelling jPas et ani201(i [Siebenmorgen et al]|20l4) we 
used other optical constants and grain size distributions. 


where 


Cext,j — 


TT/27r/27r/2 

'III 


1 F/^TM / \ I /^TE 1 

2 i^ext V^^,j 5 ■ 5 rv) J 


xfj{^.P,--)dy=>dujdP, (9) 


7r/2 TT 7r/2 

Cpoi,3i = -^J j y C'poi,si/si(C>/3, ■••) cos2il> dip duj dp , 

0 0 0 

( 10 ) 

and 

Cpoi.si = i[C7e™si(mA,si, ...,a) - Cllsi] 
for prolate spheroids and 

Cpoi.si = Ce™,si(mA,si,a) - 
for oblate spheroids. 

Here D is the distance to the star, A the wavelength, mxj, aj/bj 
and nj(rv) are the refractive index, aspect ratio and size distri¬ 
bution of spheroidal particles of the yth kind (j =Si for silicate 
particles and j =C for carbonaceous ones, respectively), rv.min.j 
and rv,max,j are the minimum and maximum radii, respectively, 
a is the angle between the wave vector of the incident radiation 
and the rotation axis of a spheroid, and the extinction 

cross-sections for two polarization modes connected with the parti- 
cle orientation relative to th e electric vector of the incident radiation 
llBohren & Huffm^ll983h . These cross-sections were calculated 
usi ng a solution to the light scat t ering p roblem for spheroids given 
by IVoshchinnikov & Farafono^ ( Il993h . The angle t/i is expressed 
through if, to, p, H (see definitions of these angles and relation s be¬ 
tween them, e.g., in lPas et al.l2010LlSiebenmorgen et al.l2014l) . and 
finally P, ...) describes the distribution of the particles of the 
jth kind over orientations. 

We calculate the extinction in the B, V bands as extinction 
A(A) (see equation|7]l averaged over wavelengths within the corre¬ 
sponding passband as follows: 


X 2 

Ax =-2.5 logio y Fx{X) 


exp 


A((A) 

1.086 


dX, 


(13) 


where Fx{X) is the normalized transmittance curve for a band X 
(=B, V) and Ai and A 2 are the band limits. 

When selecting the dust materials, we followed VH14 and 
choose the amorphous silicate with a 10% volume fraction of 
Fe (’silicate_FoFelO.RFF) and hydrogen rich aliphatic carbon (a- 
C(:H) material with a band gap Eg = 2 .5 e V). The optic al con¬ 
stants were taken from ljones et alj ( 20131) an d|I one i ( l2012h for the 
silicate and carbon, respectively. 

An important constituent of the model is the grain size distri¬ 
bution. It is obtained as a result of fitting of the observed interstel¬ 
lar extinction and polarization curves an d may be a rather compli¬ 
cated function (se elVoshchinnikovll2012l for a review). Here we in¬ 
voke t he results of lHirashitj j2012l) and Hirashita & VoshchinnikovI 
( l2014h who investigated time evolution of the grain size distribu- 
tion caused by accretion and coag ulation in an interstellar cloud. 
iHirashita & VoshchinnikovI ( 12014 ) and VH14 examined whether 
dust grains processed by these mechanisms can explain variations 
of the interstellar extinction and polarization curves observed in 
the Milky Way. They assumed that an initial grain size distribution 
(for time T — 0) should fit the mean Milky Way extinction curve 
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dWeingartner & DraindbOoH) . It was found that the observatioi 
data could be explained provided the model is ‘tuned’, i.e. wh 
coagulation of silicate dust is more efficient, with the coagulati 
threshold being removed, and coagulation of carbonaceous dust 
less efficient compared to the original model. Thus, the time 
grain processing T is uniquely related to the grain size distributit 
VH14 showed that the time scale T ~ (30 —50)(nH/10®cm“®) 
Myr, the polarization maximum shifts to longer wavelengths (An, 
grows) and the polarization curve becomes wider (K decreases), 
growth of K and Amax also occurs, when we increase the cut-( 
size rv,mt that corresponds to the interface between non-align 
and aligned grains. Both model parameters T and riz^cut influen 
the polarizing efficiency Pmax/i?(B — V) (and Pmax/^y) but t 
degree and direction of grain alignment as well as the particle ty 
and shape may produce similar effects. 

Modelling requires a specification of the function describi 
the distribution of particles over orientations according to a ; 
lected alignment mechanism. Two mechanisms are most popul 
the magnetic alignment based on the paramagnetic relaxation 
grain materi al containing about one per cent of iron impurities (I 
mechanism; IPavis & GreensteinI Il95ih . and the radiative torq 
alignment (RAT alignment) arising from an azimu thal asymmel 
of the light scattering by non-spherical particles ( IDolginov et 
I 1979 I) . Both mechanisms encounter some difficulties. 

The DG mechanism requires a stronger magnetic field th 
the average galactic one, and the polarizing grains were assum 
to contain small clusters of iron jjones & S2 itzeij[lj67[) and al 
needs to be spun up to very high velocities 1 Purcell! 197911 . 

The RAT mechanism h as been updated rath er recently (s 
Draine & WeingartneJI 19971 and the discussion in lAndersson et 


20151) . However, so far, it has been made rather at the qualitati 


level. The theory of this mechanism is not based on careful lij 
scattering calculations of interstellar polarization as the usage 
the Rayleigh reduction factor is appropriate just in the infrared p 
of spectrum (see, e.g.. lWhittet et al.ll2008h . Note also the absen 
of the correlation between the submillimeter polarization and dusi 
temperature predicted by the RAT mechanism and not de tected in 
Planck data jpianck Collaboration Xlll2014l ; lGuilletll2015h . 

Anyway, as the alignment function 
nism has not been specified in enough 
lAndersson et all l2015h , we utilized a properly modified function 
for the imperfect alignment of the spheroidal grains in Davis- 
Greenstein mechanism (IDG alignment). In this case, the distribu¬ 
tion function /3) depends on the orientation parameter ^ 

and the angle P ( iHong & Greenber j|l980h . P is the opening angle 
of the precession cone for the particle angular momentum which 
precesses around the direction of the magnetic field. The alignment 
function is written as 


for the RAT mecha- 
detail torainel |2015l; 




ismP 


(^2 cos^ p sin^ /3)3/2 


(14) 


The parameter ^ depends on the particle size rv, the imaginary part 
of the magnetic susceptibility of a dust grain x" = xWd/Td, where 
ujd is the angular velocity of the particle, hydrogen number density 
riH, magnetic field strength B, and temperatures of dust Ta and gas 
Tgas. For particles of the j-th kind. 




+ Cf (T^d.^/Tgas) 
rv^ + 


(15) 


where 


= 8.23 10^® 


niiTgasTdj 


fim 


(16) 



Figure 4. Size dependence of the normalized polarization cross-sections 
for prolate silicate spheroids with a/b = 3 (A = 0.55 fim) and the align¬ 
ment function fsp^siyP) lot P equal to its average value {/?) (T = 0, 
I'V, cut = 0.13 /rm, = 0.5 /rm). 


and we assume Tgas = 10 Td. In standard interstellar conditions, 
we have ~ 0.3 — 0 .5 /xm, which al lows one to account for the 
observed polarization llPas et al.ll2010h . Note that the DG mech¬ 
anism produces the “right” alignment when the minor grain axes 
tend to align parallel to the magnetic field if Td < Tg^s. Measure¬ 
ments show that t his condition is satisfi ed even in dense parts of 
interstellar clouds dKoumpia et alJl2oT^ . 

We suggest the following modified IDG a lignment fun ction 
with reduced alignment of small grains (see also iMatlii^l 1 9861) : 

= [l-exp(-rv/rrz,cut,j)®] x 

(17) 

where is a cut-off parameter. The alignment function 07} 

gives a smooth switch from non-aligned to aligned grains. It is 
shown in Fig. |4] for the model with the parameters: T = 0 Myr, 
fiz, cut = 0.13/rm and 5o°si = 0.5/rm (see Fig. O. The align¬ 
ment function was calculated for several grain sizes as a function 
of the mean precession cone an gle (/3) as given by equation (5) in 
lAannestad & Greenberd ( Il983h . It should be emphasized that our 
function /j (...) qualitatively agrees with that expected for the RAT 
mechanism. 

Figure |4] also shows the size dependence of the product 
Cpoi,Si(fv, «,•■•)'^Si(fv) for A = 0.55/rm normalized to its 
maximum value. We consider the non-rotating prolate spheroids 
with a/b = 3 in two orientations: a = 90° and 15°. The func¬ 
tion nsi{rv) corresponds to the initial grain size distribution for 
time T = 0. It is clearly seen that the contribution of larger par¬ 
ticles to polarization grows with decreasing the angle a between 
the wave vector and the particle symmetry axis. The figure demon¬ 
strates that a simplified treatment of the partly aligned grains using 
the Rayleigh reduction factor may be a serious error. 

Our model assumes that the magnetic field does not change its 
direction and strength in the line of sight to the star (within the re¬ 
gion where the polarization origins). However, turbulence is known 
to affect the field. A comparison of our model polarization P with 
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Polarizing efficiency and grain growth 1 


the observations may require a correction reducing the polarization 
degree as follows: Pcorrected = FP. Now one often applies the 
so called depolarization factor F that depends on the ratio of the 
regular (uniform) component Bq of the magnetic field to the total 
magnetic field Btot. = -Bo + Bt, where Bt is turbulent (random) 
component. Thus, P = 1 in the case of the regular magnetic field 
(Bt = 0) and F — 0.5 in the case of the equipartition between the 
magnetic and turbulent kinetic energy (Bt = Bq). 


4 RESULTS AND DISCUSSION 

We have performed calculations of the interstellar extinction and 
polarization curves for prolate and oblate homogeneous spheroids 
consisting of silicate and amorphous carbon. The particles of 74 
sizes in the range from r\/,min = 0.001 fim to rv,inax = 1 ftm with 
four aspect ratios a/b — 1.5, 2, 3 and 4 were utilized. The curves 
were calculated for 77 wavelengths in the range from A = 0.2 to 

5 fj,m. We computed the colour excess E{B — V), parameter Rv 
and parameters of the Serkowski curve Pmax, Amax, and K. The 
respon se functions for the B and V bands were taken from IStraizvsl 

( I 1992 I 1 . 

In this paper, we focus on a relation between the polarizing ef¬ 
ficiency Pmax/-E(B — V) and Amax. The dependence of the width 
of the polarization curve (parameter K) on the position of its max¬ 
imum has been analysed by VH14. It should be emphasized that K 
and Amax are mainly determined by the size distribution nsi(rv) 
(depending on the time of grain processing T) and the threshold 
on the size of aligned silicate grains rv.cut. At the same time, K 
and Amax are weakly affected by the degree and direction 

(fl) of the particle orientation. However, it is not the case for the 
polarizing efficiency. Therefore, we vary the alignment parameters 
and fl (keeping fixed the parameter S^q — 0.01 ^m). 

Our choice of Pmax/B(B — V) as the polarizing efficiency 
instead of Pmax/Av presents the fact that the ratio of the total 
to selective extinction Rv given by observations is not reliable 
enough. Furthermore, g rain processing in interstellar clou ds leads 
to larger values of Ry jHirashita & Voshchinnikovll2014l : see also 
Table O that produces an additional trend in the behaviour of the 
ratio Pmax/Arz (see equation[6]( with a growth of T. 

We start with a consideration of the model with prolate grains 
and a/b = 3. Variations of the grain type and shape are analysed 
in Sect. 14.41 


4.1 Dust grains without processing (T = 0 Myr) 

First, we have calculated the data for diagrams Pmax/B(B —U) vs. 
Amax for the initial size distribution of silicate and carbonaceous 
grains (Fig. [5]l. In this case, the average observational Serkowski 
curve (Amax = 0.55/im, K = 0.92) can be fitted if ry, cut = 
0.13/rm (VH14). If we assume 5™^ = 0.5 /rm, fl = 60°, and 
F = 1, the polarizing efficiency of the interstellar medium would 
achieve Pms.^/E{B — V) = 5.23 %/mag. (large filled circle in 
Fig.lDl. For F <1, this efficiency drops accordingly. 

Figure |5^ shows how variations of ry, cut influence the 
polarizing efficiency. It is clearly seen that Amax grows and 
fmax/B(B — V) decreases with a growing cut-off parameter, 
i.e. for larger values of ry, cut, the polarizing efficiency becomes 
smaller and the maximum of the curve P(A) shifts to longer wave¬ 
lengths. This effect is more pronounced for larger values of fl. 

An average size of polarizing dust grains can be found using 


the following expression: 


(t'y.poi, Si) 




/ 


rynsi(ry) dry 


^ViCut 


^Viinax 

j nsi{rv)drv 

'^V’jCUt 


( 18 ) 


For the model which gives the best approximation to the average 
observational Serkowski curve (ry, cut = 0.13/im), (ry,poi,si} = 
0.17 fim. For ry, cut increasing from 0.08 to 0.30 /tm, Amax grows, 
PraBx/E[B — V) reduces and the average grain size (ry,poi,si) 
increases from 0.12 to 0.31 pm. 

Note that for the models with ry, cut ^ 0.12 pm the polar¬ 
ization peaks in the blue part of spectrum (Amax 0.5 pm) that 
is outside the range of the observational values of Amax except for 
the case of Cyg OB2 (Fig.|5^). The models with ry, cut > 0.12 pm 
can explain only a small part of the observational data if we fix the 
parameter do^Si • Also it is evident that the cut-off parameter should 
not exceed ~ 0.2 — 0.3 pm. 

By increasing the degree of grain alignment and the angle H 
it is possible to reproduce almost all observations of open clusters 
(Fig.lSj)). However, the models do not give polarization curves with 
Amax ^0.65 pm observed in several dust clouds. In this case, even 
very large values of would not do. 

Summarizing, we can conclude that the models with the ini¬ 
tial size distribution of unprocessed dust grains fail to explain the 
observational data with large values of Amax. 


4.2 Processed dust grains (T = 10 — 40 Myr) 

Next we investigate changes of the polarizing efficiency owing to 
evolution of the grain size dist ribution. Our results are shown i n 
Fig.|6]for the tuned models (see lHirashita & Voshchinnikovll20l3) . 
The left panel illustrates the combined effects of ry, cut and H 
variations for the model with T = 20 Myr0 In this case, we 
have obtained the polarization curve that is most close to the av¬ 
erage observational Serkowski one when ry, cut = 0.095 pm, 
(ry,poi,si) = 0.15 pm {F = 1). From a comparison of Figs.|6^ 
and[^, it is clearly seen that the evolutionary effects are very im¬ 
portant: by varying ry, cut and H in the model with T = 20 Myr it 
is possible, in fact, to explain all observational data both for open 
clusters and dark clouds when B = 1, and all except for two largest 
values when F = 0.5. If we assume — 0-5 pm, the models 
with the cut-off 0.08 pm ^ ry, cut ^ 0.2 pm and H ^ 30° give 
the polarizing efficiency and Amax similar to the observed average 
values when F < 0.8. 

Note that an explanation of the available observational data 
does not require the models with more processed dust grains (T ^ 
30 Myr). This follows from Fig. [bjt where we present the models 
with increasing time of grain evolution and fixed other parameters. 
With growing T both the polarizaing efficiency and the wavelength 
of maximum polarization grow, and this effect is more pronounced 
if the direction of the magnetic field is perpendicular to the line of 
sight (i.e. for Q. = 90°). 


® This time of dust pr ocessing in a molecular cloud with the hydrogen den¬ 
sity riH = 10® cm~® jHirashita & Voshchinniko\ll2014h . 
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Figure 5. Polarizing efficiency Pmax/E{B — V) versus the wavelength of maximum polarization Amax- Symbols (filled squares and triangles and cross) 
show the average values for clouds and clusters presented in Table ^(columns 7 and 8) and Fig.[3 Vertical and horizontal dashed lines con'espond to the 
mean interstellai' value of Amax = 0.55 fim and the upper limit of the polaiizing efficiency given by equation 0. Stars (★) connected with solid and dashed 
lines show theoretical results calculated for the model: prolate spheroids, a/b = 3. Large filled circle corresponds to the model with parameters: T = 0 Myr, 
rv.cut = 0.13 Aim, = 0.5 fim, Q = 60° which is the closest to the average observational Serkowski curve (Amax = 0.55 fim, K = 0.92). Panels (a) 
and (b) illustrate the variations of the cut-off parameter ry^ cut the degree (^J^^) and direction ($1) of grain alignment, respectively. The depolarization 
factor F can be involved by a proper vertical shift of the theoretical curves. 




Amaxj P*^ Apfiax? P^^l 


Figure 6. The same as in Fig.[5]but now for the models with different time of grain evolution T. Panel (a) illustrates the variations the cut-off parameter ry^ cut 
for the model with prolate spheroids, a/b = 3, T = 20 Myr, = 0.5 /im. Large filled circle corresponds to the model with ry^ cut = 0.095 fim and 
= 45°. Panel (b) illustrates the variations of the time of grain evolution for the model with parameters: ry^ cut = 0.13 fim, = 0.5 fim, = 60°. 
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Table 2. Changes of the observed characteristics in dependence on the 
model parameters. 


Parameter 

- U) 

•^max 

K 

Ry 

cut,Si 

\ 



0 

T / 

/ 


\ 



/ 


0 

0 

n /’ 

/ 

\ 

0 

0 

prolate oblate 

/■ 



\ 

a/h (prolate) ^ 

/ 

\ 

\ 

\ 

a/b (oblate) 

/ 

0 


\ 


Note: or indicate the increase or decrease of the corresponding quan¬ 

tity with the growth of parameter. 0 means a little change of the quantity. 
Parameters: r^cut Si ~ minimum size of polarizing Si grains; T - time 
of grain processing; - degree of grain alignment; Q - direction of 

magnetic field. 


4.3 Models with prolate grains and a/b = S: summary 

Here we discuss the models in the context of simultaneous vari¬ 
ations of the four model parameters: 'rv^,cut,Si^ T, and fl. 

The first two parameters atfect the size of polarizing and processed 
grains, the last two parameters determine the degree and direction 
of grain alignment. The type and shape of the particles (prolate 
spheroids with a/6 = 3) remain the same in all the models. 

If we assume that the polarizing efficiency on average is 
Pins.x/E{B — V) =4—6%, the average observational Serkowski 
curve can be reproduced by the models with Q = 60° (T = 0 Myr) 
and = 45° (T = 20 Myr) for F = 1. These models are plotted 
by large filled circles in the left and right panels of Fig.|7] respec¬ 
tively. Using the circles as starting points, we have calculated a set 
of models with reasonable variations of the four mentioned param¬ 
eters. The final picture thus obtained resembles an octopus shown 
in Fig. [7] where just one model parameter varies along each line. 
A distinguishing feature of the Figure is the different behaviour of 
these four lines, which allows one to estimate the effects of changes 
of the parameters or their combination. Evidently, the additional in¬ 
formation about the width of the polarization curve (the parameter 
K of the Serkowski curve) would be helpful. Variations of the ob¬ 
served polarization with changes of the model parameters show the 
trends summarized in Table For F < 1, the octopus moves to 
the bottom accordingly. 

Figure[7^ demonstrates that the models with unprocessed dust 
grains (T = 0 Myr) cannot explain the observational points in the 
right and upper parts of the plots Pms.^/E{B — V) vs Amax- The 
models with a higher value of the alignment parameter <5™^ do not 
solve the problem. Thus, we need to consider the models with larger 
dust grains resulting from dust evolution in molecular clouds (see 
Fig-0’)- As follows from this figure, the models with T = 2Q Myr 
permit to reproduce the average observational data for clouds and 
clusters under consideration if we vary three additional parameters: 
fv,cut,Si, fl. This fact illustrate the models plotted by 

dot-dashed lines in Fig. [TJj for Sjfgt = 1.0 and two values 
of the depolarization factor F. The upper curve shows the results 
in the case of the regular magnetic field (F = 1) while the lower 
curve is obtained for the model that takes into account the turbulent 
component of the magnetic field (F = 0.5). 

A detailed analysis of observational data for separate objects 
must be made more carefully involving information on interstellar 
extinction and will be carried out in the next paper. 


4.4 Variations of grain type and shape 

The interstellar polarization is often stated to p rovide informa tion 
about the shape of cosmic dust grains (see, e.e.. lWhittetll2003l) . So 
far such modellin g has been made for infinite circular cylinders 
and spheroids (see lVoshchinnikovll2012l for a review). The infinite 
cylinders are the simplest non-spherical particles whose shape can¬ 
not be varied. The spheroids were con sider ed just in a few papers . 
In particular, [Rogers & Martini (Il979h and iKim & Martini ll 1995ll 
used prolate and oblate spheroidal grains with different aspect ra¬ 
tios a/6. But their consideration was restricted by perfectly aligned 
rotating/non-rotating particles with a fixed orientation of the mag¬ 
netic field relative to the line of sight 90°. However, in this case 
the polarizing efficiency should be in several times larger than the 
observed maximum even for F = 0.5. 

We have calculated the interstellar extinction and polariza¬ 
tion for prolate and oblate spheroids with four aspect ratios a/6 = 
1.5, 2, 3 and 4. Some of these results are given in Fig.[^and Table[3] 
for rv.cut = 0.1, 0.2 pm and T = 0, 20 Myr. Comparing Fig. [8] 
with Fig.[7]we can arrive to the conclusion that for prolate spheroids 
variations of the aspect ratio seem to be less important than varia¬ 
tions of other model parameters, namely rw.cut ,Si, T, 6™^ and fl. 
When the parameter a/6 grows, Amax reduces for prolate spheroids 
but remains almost the same for oblate spheroids. Note also that the 
oblate particles are more efficient polarizers and can easily explain 
extremely large values of the polarizing efficiency observed, for ex¬ 
ample, in Musca (Fig.0) when F ~ 0.8. 

A difference between prolate and oblate spheroids also ap¬ 
pears in the shape of the polarization curve: a growth of the param¬ 
eter a/6 makes the curve wider (K reduces) for prolate prides 
and narrower {K increases) for oblate particles (Table[3[j VH14 
compiled the parameters of Serkowski curve for a sample of 160 
lines of sight and found the observational limits: K ~ 0.5 — 1.5 
and Amax ~ 0.35 — 0.8 pm. As follows from Tablej^ the major 
part of the models with oblate particles is beyond these limits. 

A more detailed consideration of the shape effects is outside 
the scope of this work and requires a further investigation. 

4.5 General discussion 

It is evident that any modelling allows determining just of in¬ 
tervals for the model parameter values. Nevertheless, it provides 
some trends in the behaviour of the observed characteristics (Ta¬ 
ble [2}. Using the diagrams Pmax/i5(f3 — V) vs. Amax (Figs. |5] 
- [8} and K vs. Amax (VH14) we can conclude that our model 
has two key parameters: the threshold on the size of aligned sil¬ 
icate grains ri/,cut,Si and time of grain processing T. After com¬ 
paring with observations, it is possible to restrict the parameters: 
0.05 pm rw,cut,Si ^ 0.2 pm and T 20 — 30 MyrQ Note 
that a growth of T is accom panied by an increase of the total t o se¬ 
lective extinction ratio Rv jHirashita & Voshchinnikovll2014l : Ta- 
ble[3. Variations of the key parameters would explain the general 
observational trends while additional polarization details could be 
reproduced if we changed degree and direction of grain alignment 
6™^ and Q.. Calculations show that we can hardly distinguish the 

These tendencies are noticed in Table |2] 

^ This time does not contradict to the lifetime of molecular clouds 
jPagani et alJl201l1 : lDobbs & Pringle ||2013l) because the time-scales of ac¬ 
cretion and coagulation oc jHirashita & Voshchinniko\j|2014l) . i.e. if 
we adopt njj = 10'* cm“^ instead of 10® cm“®, the same size distribu¬ 
tions are reached at one-tenth the time. 
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^max5 ^maxs 

Figure 7. Polarizing efficiency Pmayi/E{B — V) versus the wavelength of maximum polarization Amax- Symbols show the observational data presented in 
Table ^(columns 7 and 8) and Fig. [5] Stars (★) connected with solid and dashed lines illustrate theoretical results calculated for the model: prolate spheroids, 
a/b = 3, T = 0 Myr (a) and T = 20 Myr (b). Large filled circles correspond to the models closest to the average observational Serkowski curve (see Fig.[5)- 
Circles connected with dot-dashed lines illustrate theoretical results calculated for the models with T = 20 Myr and the depolarization factors F = 1 (filled 
circles) and F = 0.5 (open circles). Panels illustrate the variations of the different model parameters. 



Amaxj 



Amax? 


Figure 8. Polarizing efficiency Pmax/E{B — V) versus the wavelength of maximum polarization Amax- Symbols show the observational data presented in 
Table^(columns 7 and 8) and Fig. [3] Stars (★) connected with solid and dashed lines demonstrate theoretical results calculated for the model: = 0.5 fim, 

= 60°, T = 0 Myr (a) and T = 20 Myr (b). Panels illustrate the variations of the particle type (prolate/oblate) and particle shape (aspect ratio a/b). The 
values of the parameters ry, cut ^^d a/b are indicated. 
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Table 3. Observed chai'acteristics in dependence on the particle shape and time of grain processing T (the model with <5g®gJ = 0.5 ^m, O = 60°). 




T = 0 Myr 


T 

= 20 Myr 



ajh 

Fmax/ECB 

— ^max 

K Rv 

Pmax/£(S - V) 

'^max 

K 

Rv 


Prolate 

spheroid, rv,cut 

= 0.1 //m 

Prolate spheroid, ry^cMt 

= 0.1 fj.m 


1.5 

2.96 

0.518 

1.135 3.633 

4.46 

0.612 

1.013 

6.060 

2.0 

4.88 

0.501 

1.073 3.580 

7.23 

0.595 

0.969 

5.971 

3.0 

7.10 

0.472 

0.983 3.397 

9.91 

0.562 

0.904 

5.534 

4.0 

8.44 

0.446 

0.925 3.218 

10.84 

0.533 

0.861 

5.045 


Prolate 

spheroid, rv',cut 

= 0.2 //m 

Prolate spheroid, rv,cut 

= 0.2 fim 


1.5 

1.21 

0.689 

1.370 

2.956 

0.843 

1.326 


2.0 

1.93 

0.667 

1.282 

4.506 

0.816 

1.246 


3.0 

2.68 

0.627 

1.161 

5.708 

0.770 

1.136 


4.0 

3.05 

0.594 

1.088 

5.947 

0.729 

1.070 



Oblate 

spheroid, ry^cut 

= 0.1 

Oblate spheroid, rv,cut 

= 0.1 fim 


1.5 

5.62 

0.567 

1.357 3.495 

8.18 

0.679 

1.050 

5.708 

2.0 

9.23 

0.580 

1.442 3.384 

12.29 

0.697 

1.107 

5.146 

3.0 

13.72 

0.584 

1.524 3.197 

15.73 

0.703 

1.194 

4.266 

4.0 

16.95 

0.580 

1.681 3.109 

17.72 

0.695 

1.306 

3.822 


Oblate 

spheroid, ry^cut 

= 0.2 

Oblate spheroid, rv,cut 

= 0.2 fim 


1.5 

2.73 

0.763 

1.700 

7.05 

0.983 

1.297 


2.0 

4.87 

0.784 

1.877 

11.93 

1.004 

1.489 


3.0 

8.02 

0.794 

2.035 

17.01 

0.996 

1.820 


4.0 

9.91 

0.776 

2.242 

19.07 

0.969 

1.992 



cases when the direction of the magnetic field is near the line of 
sight (O ^ 15°) or close to perpendicular to the line of sight 
(fl = 75° — 90°). A distiction of the degree of grain alignment 
also presents a problem if ^3 — 5 pm. 

The overall picture of variations of the polarizing efficiency 
presented above would not dramatically change if we varied the as¬ 
pect ratio a/b of prolate grains (Fig. [8] Table[^. Thus determining 
the shape of the particles like ‘cigars’ would be problematic. 

The polarizing efficiency for oblate particles significantly ex¬ 
ceeds that of prolate ones. Thus, it is possible to reduce the align¬ 
ment degree or particle aspect ratio. However, very flattened ‘pan¬ 
cakes’ (a/b ^ 3) must be ruled out because they produce too nar¬ 
row polarization curves. 


5 CONCLUSIONS 

The main results of the paper can be formulated as follows: 

(i) We collected the observational data about the polarizing effi¬ 
ciency Pmax/E{B — V) and the position of the maximum polar¬ 
ization Amax for 243 stars in 17 objects separated into two groups: 
dark clouds and open clusters. The average data for these two 
groups are distinguished by the parameter Amax. For open clusters, 
the mean values of Amax are grouped around the average interstel¬ 
lar value 0.55/rm with rather small deviations. For dark clouds, the 
values of (Amax) may be significantly larger than 0.55/rm and have 
a wide scattering. 

(ii) We utilized th e mo del applied earlier by 
IVoshchinnikov & Hirashital ( 12014|) and including time evo¬ 
lution of the grain size distribution due to the accretion and 
coagulation processes and new optical constants of grain materials. 
To calculate the polarization curves, we used homogeneous silicate 
and carbonaceous spheroidal particles of different aspect ratios a/b 
having imperfect alignment. It was assumed that polarization was 
mainly produced by large silicate particles with sizes ry ^ r^^cut. 


We calculated the wavelength dependence of extinction and 
polarization and determined the parameters PmBx/E{B — V) and 
Amax. 

(iii) We focused on the interpretation of the observed relation 
between Pinax/E{B — V) and Amax keeping in mind probable de¬ 
polarization effect caused by the random component of the mag¬ 
netic field. Theoretically both quantities are mainly determined 
by two key parameters: the threshold on the size of aligned sili¬ 
cate grains rw,cut and time of grain processing T. We found that 
the models with the initial size distribution (no grain processing, 
T = 0 Myr) that reproduce the average curve of the interstellar 
extinction fail to explain the data with Amax ^0.65 pm observed 
in several dark clouds. An inclusion of evolutionary effects allows 
us to find the models which give the polarizing efficiency and Amax 
values similar to the average observed ones for reasonable depolar¬ 
ization (F ~ 0.5 — 0.8). This occurs for 0.05 pm rv',cut,Si 
0.2 pm and T (20 — 30)(nH/10®cm“®)“^ Myr and is in agree¬ 
ment with the lifetime of molecular clouds. The additional polar¬ 
ization details are reproduced if we change the degree (5o®®) and 
direction (H) of particle orientation. However, these results can¬ 
not yet be a conclusive evidence that the dust evolution only can 
completely explain the observational data under consideration be¬ 
cause of a rather large number of the parameters still involved in 
our modeling. A more detailed analysis in particular that involving 
the width of the polarization curve is required. 

(iv) We found that generally the change of the aspect ratio a/b 
of prolate grains would not dramatically influence the variations of 
the polarizing efficiency that makes the determination of the grain 
shape rather problematic for prolate particles. The polarizing ef¬ 
ficiency for oblate spheroids significantly exceeds that of prolate 
ones but particles with a/b ^ 3 must be excluded because they 
produce too narrow polarization curves that are not observed. 
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